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I. INTRODUCTION 

Quantum electrodynamics (QED) in canonical form 
can be formulated in a number of different ways, because 
the QED Lagrangian is independent of the velocity of 
the scalar potential. Formulations can be broadly cat- 
egorized into two types. The first type comprises the 
manifestly covariant formulations within a Lorcnz-type 
gauge such as the Feynman gauge. These require the 
addition of a gauge symmetry-breaking term to the La- 
grangian, after which the scalar potential can be taken 
as a legitimate canonical variable with a nonvanishing 
conjugate momentum. Quantization is carried out using 
an indefinite metric with the Lorenz subsidiary condition 
implemented as a constraint defining the physical space 
of states [l[ . 

The second type of formulation consists of those in 
which the scalar potential is eliminated from the La- 
grangian altogether. The vector potential is then the 
sole canonical coordinate in the Hamiltonian. Of course, 
the relative simplicity of such an approach comes at the 
price of sacrificing manifest covariance Q. 

In the Weyl gauge the scalar potential is set identi- 
cally equal to zero. In the Coulomb and Poincare gauges 
it is eliminated in favor of the gauge-invariant degrees 
of freedom. How this is achieved, varies between the 
two gauges, so the resultant Lagrangians are distinct. In 
these cases the gauge is fixed at the classical level with 
a particular Lagrangian belonging to a particular gauge 

ai. 

Conventional calculations in atomic physics and quan- 
tum optics most often start with a nonrelativistic Hamil- 
tonian in the Coulomb or Poincare gauge With 
regards to such calculations there are two main draw- 
backs. The first is their nonrelativistic description of the 
material degrees of freedom. The second is their ab ini- 
tio restriction to a particular gauge, which makes them 
incapable of fully exploring the consequences of gauge 
freedom. Indicative of the potential benefit of a more 
general formulation is the fact that certain admixtures 
of the Coulomb and Poincare gauge Hamiltonians have 
proven useful in quantum optics 



In this paper we focus on formulations of QED of the 
second (noncovariant) type in which the scalar poten- 
tial is eliminated. Specifically the aim of this paper is 
to address both of the issues raised above; first using a 
quantized Dirac field description for the material degrees 
of freedom and second keeping at the forefront the gauge 
freedom present in the theory. 

The latter is achieved following Ref.[l(| by reexpress- 
ing the gauge-dependent longitudinal degrees of freedom 
using the gradient of a linear functional of the gauge- 
invariant transverse degrees of freedom. The functional 
is defined in terms of an arbitrary c-number function 
g(x, x'), which then carries the gauge freedom of the the- 
ory. 

As a description of the Maxwell field we use a Hilbert 
space of Schrodinger wave functionals of the vector po- 
tential in which an inner product can be defined through 
functional integration [ll | . The material degrees of free- 
dom can be defined similarly in terms of a Hilbert space 
of functionals of a Grassman field [ll| . 

To obtain a Hamiltonian in an arbitrary noncovariant 
gauge (defined by g) we combine these elements with 
an adaptation of the quantum- mechanical gauge-fixing 
method presented in [12|]- In this approach one initially 
adopts the Weyl gauge. From there the physical subspace 
of states (wave functionals) is defined as the subspace of 
states vanishing under the action of the Gauss law con- 
straint. The Hamiltonian can be found in different non- 
covariant gauges through the use of unitary gauge-fixing 
transformations acting on the Weyl gauge Hamiltonian. 

Particular choices of the function g can subsequently 
be made to render the Hamiltonian in a fixed gauge. Two 
particular choices yield the Coulomb gauge and Poincare 
gauge Hamiltonians. 

On the one hand our results are of importance for 
the relativistic theory of atoms and molecules and on 
the other they provide an interesting extension of the 
quantum-mechanical gauge-fixing formalism already em- 
ployed in noncovariant QED. We will show that this ex- 
tension has important implications with regards to the 
ontology of noncovariant QED in general, because it al- 
lows for a thorough exploration of the gauge freedom of 
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the theory in canonical (Hamiltonian) form. The require- 
ment of gauge invariance of a result is translated into the 
requirement that it be independent of the choice of the 
arbitrary function g. 

There are four sections to this paper. In Sec. [H]we use 
the quantum-mechanical gauge-fixing formalism to ob- 
tain a Hamiltonian in an arbitrary noncovariant gauge. 
In Sec. Mil we discuss the implications of the gauge free- 
dom still present in our formulation. We then define 
useful gauge-invariant operators and as an application 
address the problem of causality in spatially separated 
material systems [l3|. In Sec. IIVI we finish with a brief 
summary of our results. 



II. THE QED HAMILTONIAN IN AN 
ARBITRARY NONCOVARIANT GAUGE 



We will begin in Sec. Ill Al bv outlining how the general 
idea of gauge fixing we are going to employ works in the 
simple case of classical electrodynamics. The aim is then 
to extend this idea to noncovariant relativistic QED. To 
this end we start Sec. Ill Bl formally with the QED La- 
grangian and obtain from it the Hamiltonian and canon- 
ical operators, which are supposed to satisfy canonical 
commutation relations. We go on to identify the Gauss 
law constraint, which defines the physical subspace of 
states and review a particular class of transformations 
called residual gauge transformations [l2[ • 

In Sec. Ill CI we identify the states of the system as 
Schrodinger wave functionals and determine the general 
form of a physical state using the coordinate represen- 
tation for the canonical operators of the Maxwell field 
and the Gauss law constraint. From there we identify a 
general unitary gauge-fixing transformation U s as a map 
from the physical space of states H p to a space % gl which 
is the space of states for the gauge g. Next in Sec. Ill Dl 
we determine the effect of this transformation on the var- 
ious operators of the theory and express the Hamiltonian 
in the arbitrary gauge g. Finally in Sec. Ill El we use the 
Hamiltonian to calculate the Dirac equation in the gauge 



A. Gauge in classical electrodynamics 

In electrodynamics the electric and magnetic fields de- 
fined by 



E = -V0- 



dA 

~dt' 



B = V x A (1) 
are invariant under a change of gauge; 



A' = A + V/, 

= +-8i 



(2) 



for an arbitrary function /. Clearly the transverse vec- 
tor potential At is gauge- invariant, while the (redun- 
dant) gauge- dependent degrees of freedom are longitudi- 
nal. For this reason the transverse vector potential serves 
as a convenient starting point relative to which vector po- 
tentials in other gauges can be defined by 



V/, 



(3) 



with the function / determining the gauge. The gauge- 
invariant transverse electric field is given by 



Ei 



while the Gauss law 



VE = p 



(4) 



(5) 



involving the charge density p ensures that the longi- 
tudinal electric field is equal to minus the gradient of 
the (static) Coulomb potential; El = —W. Using this 
equality and Eq. ([1]) we see that the scalar potential 
accompanying the vector potential in Eq. ([3]) can be 
written 



at 



(6) 



The Coulomb gauge is defined by the choice / = 0, but 
it is not the only way in which the longitudinal degrees of 
freedom can be eliminated. In nonrelativistic QED it has 
been shown [lj| that a general gauge-fixing condition is 
given by a linear functional constraint satisfied by vector 
potentials of the form 



A = A T + V J dVg(x, x') • At(x') 



(7) 



where g(x, x') is the Green's function for the divergence 
operator; 



V-g(x,x') = <5(x-x') 



(8) 



In Eq. ([7]) the redundant degrees of freedom have been 
re-expressed through a functional of the gauge-invariant 
transverse degrees of freedom. 

While the longitudinal component of the Green's func- 
tion g is fixed according to Eq. © by 



g L (x,x') = -V 



4ir x — x' 



(9) 



its transverse component is essentially arbitrary, meaning 
that the gauge is determined through a choice of gT- This 
idea has been employed in nonrelativistic QED to obtain 
a Hamiltonian in an arbitrary gauge [lOj , but it has yet to 
be extended to the relativistic setting. Furthermore, the 
nature of gauge transformations and residual symmetries 
in such a framework has not been explored. 



3 



B. The QED Lagrangian and Hamiltonian 



We start formally with the QED Lagrangian density 

1 



if = i-fol^D^j - ( 7o m + e<j> e )x^xp - -^F^ 



(10) 



where D„ = 9„ + ie.A M is the gauge covariant derivative, 
= dfj,A v — d v A p is the electromagnetic field strength 
tensor and <p e is an external potential due, for example, 
to nuclei. Since the Lagrangian is independent of the 
velocity of the scalar potential its conjugate momentum 
is identically zero. As a result it is natural to quantize the 
theory within the Weyl gauge corresponding to the choice 
(j) = 0. The remaining redundant degrees of freedom are 
eliminated by defining the physical subspace of states H p 
consisting of those states, which vanish under the action 
of the Gauss law constraint G \ip p ) = (V • E — p) \ip p ) = 0. 

The Hamiltonian density is obtained from the La- 
grangian density via a Legendre transformation 

,34? = -i^ot ■ (V - ieA)tp + ((3m + e(j> e )^ip 

+ i[n 2 + (VxA) 2 ]. (11) 

quantum-mcchanically tp and its conjugate are Dirac 
field operators satisfying the anti-commutation relation 



{^(x),V t (x')} = 5(x 



(12) 



while A and II = — E are the canonical operators of the 
Maxwell field satisfying the commutation relation 



[A i (x),Il j (x')]=i5 ij 6(x-x') 



(13) 



We employ the usual definitions of the charge and cur- 
rent densities 



p = eip'ip, 
j = etp^atp, 



(14) 



in terms of which the (conserved) Noether four-current 
is j 1 * = (p, j). Related to the charge density is the polar- 
ization field P g defined by 



P g (x) = - / dVg(x,x>(x'), 



whose longitudinal component satisfies 
-V-P g = p, 



(15) 



(16) 



but whose transverse component is arbitrarily deter- 
mined by gT- Equation (|16l) bares a close resemblance 
to the Gauss law G, which can indeed be written 



g = v-u + P = v-(n-p g ). 



(17) 



As an operator G is a symmetry of the Hamilto- 
nian [G, H] — and is responsible for generating time- 
independent gauge transformations of the vector po- 
tential and Dirac field operators. Identifying a group 



({/3(x)}, +) consisting of real- valued functions on R 3 and 
group operation of addition, we define a group action $ 
acting on the vector potential and Dirac field by 



$[A,p] = A + V/3. 



(18) 



The action is implemented through unitary transforma- 
tions ft generated by g, viz., iltpft^ 1 = $>[ijj,0\ and 
QAft- 1 = <5>[A,P], where 



fif 



exp ( i / d 3 x (II ■ V + /o)/?(x) 



(19) 



These transformations are called residual gauge transfor- 
mations, with the word residual intended to signify that 
the above time independent symmetry is what remains 
of the loc al g auge symmetry present in the original for- 
mulation fial. 



C. Unitary gauge-fixing transformations 

We now turn our attention to the procedure of gauge 
fixing. In order to determine the form of a general gauge- 
fixing transformation we first need to identify the form 
of a physical state. To do this we take as a Hilbert space 
% for the composite system wave functionals tp [A] of the 
c-number vector potential A, which take values in the 
Hilbert space of the Dirac field operators fllj ]. 

A realization of the algebra of the Maxwell field oper- 
ators A and II is given on H using the coordinate repre- 
sentation 



(A^)[A] = A^[A] 



(ftp) [A] = - 



MM 



(20) 



where we have introduced carets to distinguish between 
operators and c-numbcr vector fields. Defining a scalar 
function a by Va = Al , we can vary the wave functional 
tp with respect to a and make use of Eq. ([20)1 to obtain 



.Sip 
5a 



-V ■ Iltp. 



(21) 



Using the constraint G in Eq. (jTT)) . we get for a physical 
state ip p 



= ptp p 

oa 



(22) 



and finally solving this equation gives the general form 
of a physical state; 

tp P [A] = tp p [A T + Va] 

= exp^-z J d 3 xa{x)p(x)^(pp{A T } . (23) 

Having determined the form of a physical state, we 
can begin to define some unitary gauge- fi xing transfor- 
mations. In the original work of Lenz et al. [12( a unitary 
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gauge-fixing transformation yielding the Coulomb gauge 
representation was given as 



^exp^/d^xMx)) 



(24) 



where d is defined analogously to a by Vd = Al. In 
the present context we see clearly that U eliminates the 
dependence of the physical state on Al; 



[U<pp)[A] = V P [A T ] 



(25) 



Now, as in Eq. (0 of Sec. Ill Al we write the longitu- 
dinal vector potential as the gradient of a functional of 
the transverse vector potential; 



A L = V Xg (x, [A T ]) 



(26) 



where 



Xg (x, [A T ]) = J dVg(x,x') ■ At(x') . (27) 

We can then define a more general unitary gauge-fixing 
transformation U g by 

U g ee exp (i J d 3 x |d(x) - Xg (x, [A T ])}p(x)J , (28) 

mapping from H p to a space denoted T-L g , which is the 
space of states for the gauge g; 

(U K <p p )[A] = exp^-i J d 3 xx g (x, [A T ]) ( o(x)^ p [A T ] 
= ip p [A T + Vxg] = <y5 g [A T ] G H g . 

(29) 

The (transverse component of the) Green's function g 
is essentially arbitrary and determines the gauge. Two 
commonly used examples are the Coulomb gauge gT = 
and the Poincare gauge gT(x, x') = — f Q dAx'5(x — Ax') 
0, H| • The vector potential operator in the gauge g is 
A(x) = At(x) + Vx g (x, [At]) with action on H g given 
by 



(A^ g )[A T ] = (A T + VxgVg[A T ] 



(30) 



Finally we define a unitary transformation from a fixed 
gauge g to a fixed gauge g' by 

U ss > = 

exp^i J d 3 a;{xg(x, [A T ]) - X g ' (x, [ a t])}p(x)^), 

(31) 

an example of which is the well known Power- Zienau- 
Woolley transformation 0fl] used to obtain the Hamil- 
tonian in the Poincare gauge from the Hamiltonian in the 
Coulomb gauge. Such a gauge transformation is not to 
be confused with the residual gauge (symmetry) trans- 
formation given in Eq. (|19[) . 



D. The Hamiltonian in the gauge g 

To obtain the Hamiltonian in the gauge g we need to 
determine the effect of the transformation in Eq. (|28p 
on the various operators of the theory, namely, ip, A 
and II. In doing so we will resume denoting operators 
without carets. Clearly U g leaves the vector potential A 
unchanged, while the action of A on 7i g is given in Eq. 
(J3UJ) . The effect on the Dirac field operator ip is that of 
a gauge transformation 



u K i>u- 



Je(a-x s ) 



(32) 



The canonical momentum II transforms as 

r-l TT T3 (33) 



u g uu-' = n + p g , 



so that in the new representation II represents the (neg- 
ative of) the gauge-dependent displacement operator 
D g = E + P g . Using Eq. (fTBf we find the constraint 
G and the residual gauge transformation f2 transform as 



umu: 



u g n[(i]uz 



v n, 



exp (i J d 3 x (n ■ V)/3(x)^ , 



(34) 



which arc independent of the gauge g. The constraint G 
implies that the longitudinal canonical momentum IIl 
vanishes on % g . On the one hand this means Pl alone 
represents (the negative of) the longitudinal electric field 
and on the other that the Hamiltonian density on H g can 
be written in terms of the transverse operators At and 
n T only; 

,ye = -iifta ■ [V - ie(A T + Vx g )] i' + ifim + e0 e )V>V 

+ ^(Pl) 2 + \ [(H T + P|) 2 + (V x At) 2 ] (35) 

where At + Vx g is simply the vector potential A in 
the gauge g. Equation ([33]) gives a Hamiltonian in an 
arbitrary gauge, which is fully relativistic in the material 
degrees of freedom. It is one of the main results of this 
paper. 

The commutator of the transverse operators follows 
from Eq. (|13[) and is given by 



[A T , i (x) ) n TJ (x')] = ^(x-x'), 



(36) 



with S T denoting the transverse delta function. We note 
also that by denoting the Fourier transforms of At and 
IIt with tildes we can define photon creation and anni- 
hilation operators in the usual way 



ax(k) 



wI T ,A(k) +in T , A (k) 



(37) 



where A = 1, 2 denotes one of two polarization directions 
orthogonal to k. The bosonic commutator 



[a A (k),4,(kO] = <J AA ,<5(k-k'), 
follows from Eq. (|36l) . 



(38) 
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E. The Dirac equation in the gauge g 

It is an instructive exercise to calculate in the arbitrary 
gauge g the equation of motion for the Dirac field opera- 
tor ip, which should be the Dirac equation in the presence 
of a Maxwell field. The calculation demonstrates how the 
scalar potential, like the longitudinal vector potential, is 
rccxprcsscd through the functional Xg- 

Writing first the products of Dirac field operators ap- 
pearing in Eq. (|35[) in normal order, we obtain 



zip = ja- [ - iV - e(A T + Vxg)] +/3m + e^ e 
P(x') 



e 

47T 



dV 



X - x' 



i y'dVgT(x ) x')- [n T + pf]J 



V> ■ (39) 



The first term on the second line of Eq. (f3"9")l is equal 
to eV, with V denoting the static Coulomb potential of 
charges. The transverse electric field in the gauge g is 
Et = — (IIt + P T ) an d it is straightforward to verify 
that Ex = -At, as in Eq. These equalities imply 
that Eq. (|39f is indeed the Dirac equation 



iip = [a • ( - iV - eA) + /3m + e(0 e + 4>)] ip (40) 
where we have defined the scalar potential anew by 



4> = v 

which is analogous to Eq. 



dt ' 



(41) 



III. SOME IMPLICATIONS OF THE 
FORMALISM 

Having obtained the Hamiltonian in an arbitrary 
gauge, we discuss in this section some implications re- 
sulting from the freedom to choose gx- In Sec. IIII Al we 
point out that as in the nonrclativistic case a canonical 
partitioning of the Hamiltonian is gauge-dependent and 
therefore leads to gauge-dependent definitions of quan- 
tum subsystems [lfj |. 

By means of analogy with classical electrodynamics, in 
Sec. IIII Bl we demonstrate how gauge- invariant subsys- 
tem components might be defined. Finally in Sec. IIII CI 
we consider an application of such a definition in the con- 
text of energy transfer and causality in Fermi's two-atom 
problem. 



A. Ambiguity in defining quantum systems 

To understand what the arbitrariness of gx might 
mean with regard to the physical predictions of our for- 
mulation it is important to identify the physical observ- 
ables in a given gauge g. In the Weyl gauge the canonical 



momentum II is the negative of the electric field E. In 
the gauge g we have E = -U^UU^ 1 = -(II + P g ). The 
operator IIx + P T appearing in Eq. (|3"5j) therefore rep- 
resents the negative of the transverse electric field Ex- 
This identity was used in obtaining the Dirac equation 
CH). 

Due to the gauge dependence of P T , the operator IIx 
is implicitly gauge-dependent in that it represents a dif- 
ferent physical observable in each different gauge. Ex- 
plicitly rix represents the transverse component of (the 
negative of) the gauge-dependent displacement operator 
D g . 

Now, the most common way to use a Hamiltonian (den- 
sity) such as Eq. (|35|) is to split it into free and interacting 
components as follows; 



where 



■j4?d 



3Cd + ^?M, 

—iip 1 * a ■ V"0 + (/3m + e^^ip^tp 

\2 , 1 /Tjg\2 



(42) 



i 



M 



i(P L ) 2 + i(P|)S 
[n T + (V x A T ) 2 ] 



and 



= -e^a.- (A T + Vxg)^ + IIx-Pi 



(43) 



(44) 



The reason for this splitting is of course that the sets of 
operators and {Ax, IIx} are mutually commut- 

ing. The component J%h represents the Dirac field sub- 
system, J4?m represents the Maxwell field subsystem, and 
Ji?i represents their interaction. The problem with such 
a splitting is that the subsystem components so defined 
are like the operator IIx, implicitly gauge- dependent and 
as a result physically ambiguous. 

Given this ambiguity it is natural to try and determine 
what kind of calculations can be carried out that yield re- 
sults independent of g, i.e., that are gauge-invariant. For 
example, S'-matrix elements on-energy-shell are gauge- 
invariant to all orders in perturbation theory [1J, Ha ]. 
Consequently, the conceptual difficulty regarding the def- 
inition of subsystems does not effect scattering theory 
in any way. The invariance of S'-matrix elements rests 
on the fact that the bare states (eigenstates of the free- 
energy operator Hq) coincide asymptotically with eigen- 
states of the total Hamiltonian, that is, that bare states 
are asymptotically stable. 

In order that conventional (perturbativc) calcula- 
tions using bare states produce gauge-invariant results 
a condition of free-energy conservation must be satis- 
fied [l4[. Otherwise calculations will in general yield 
gauge- dependent results. Free-energy conservation is a 
condition that must be imposed from outside the the- 
ory, so various approximations, which essentially ensure 
free-energy conservation by giving rise to a (5-function in 
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the initial energy minus the final energy of the process 
under study, are used throughout nonrelativistic QED 
and quantum optics; they include the resonant state on- 
energy-shell approximation fLU [l6j . the pole approxi- 
mation [13], the Markovian approximation (l7j . and the 
Fermi approximations 0, l| ■ 

A general survey of the types of calculation used in 
practice in nonrelativistic QED and their dependence 
on g can be found in [l4|- An example of the differ- 
ence in predicted results from calculations in two differ- 
ent gauges is well known in nonrelativistic QED where 
the Coulomb gauge (minimal coupling) Hamiltonian and 
Poincarc gauge (multipolar) Hamiltonian yield different 
results for, among other things, the theoretical lineshape 
of spontaneous emission [H . The source of this difference 
lies in the use of physically different canonical operators 
in determining the lineshape observable. 

The ambiguity regarding subsystems defined using 
canonical variables is also present at the classical level. 
It occurs when moving to the Hamiltonian formalism 
from the Lagrangian formalism. Since canonical mo- 
menta arc defined in terms of the Lagrangian, equivalent 
Lagrangians yielding equivalent Hamiltonians will in gen- 
eral not yield physically equivalent canonical momenta. 
Again, a well known example is given by the Coulomb 
gauge and Poincare gauge formulations of classical elec- 
trodynamics. 

At the quantum level the ambiguity in the defini- 
tion of subsystems can be viewed as a generic trait of 
interacting theories, whether they are relativistic, field 
theoretic, or otherwise. Given a Hamiltonian depen- 
dent on two sets of mutually commuting operators {x{\ 
and {yi} and a splitting of the Hamiltonian of the form 
H = H x ({xi}) + H y ({yi}) + H xy ({xi},{yi}), an equiv- 
alent Hamiltonian is obtained by a unitary transforma- 
tion H' = UHU -1 . In general, the subsystem compo- 
nents of H' will not be equivalent to those of H, i.e., 
H' x ^ UHxlI" 1 , with the same being true for H y and 
H xy . The importance of this fact for concepts such as 
quantum entanglement and decoherence has been recog- 
nized in the philosophy literature [l8j . 



B. "Classical- type" gauge-invariant subsystems 
and operators 

In the preceding section we discussed the gauge depen- 
dence of splitting the Hamiltonian into free and interac- 
tion components and we reviewed the situation regarding 
the kinds of calculation that yield gauge- invariant results. 

In this section we address the complimentary question 
as to whether manifestly gauge-invariant subsystem com- 
ponents can be defined from the outset. This question 
has so far received no direct attention in the literature, 
but it is important if one wishes to identify how gauge- 
invariant results might be obtained outside the range of 
validity of the simplifying assumptions of scattering the- 
ory. 



To see how we might define gauge-invariant subsys- 
tems consider first the case of a free classical electron 
coupled to a classical Maxwell field. For this system the 
Hamiltonian in any gauge can be written 

H = l -mv 2 + X - j d 3 x [E(x) 2 + B(x) 2 ] , (45) 

which represents the total energy of the system as the 
sum of the kinetic energy of the electron and the energy of 
the electromagnetic (EM) field. Regarding the electron 
variables, what varies between gauges is the identifica- 
tion of the gauge-invariant velocity r, with the electron 
canonical momentum p = mr — eA(r) 0, Holll9j|. 

The classical electron velocity and canonical momen- 
tum have clear analogs in QED; the velocity density 

f = -iip f a- [V-ie(A T + Vxg)]^ (46) 

is manifestly gauge-invariant, while the canonical mo- 
mentum density = —iifj^a ■ V?/> depends on g. Analo- 
gously to Eq. (|4"5"|) we can split the Hamiltonian into two 
gauge-invariant components 

Mr — Mm + mem, 
m° m = -iip^a- [V-ie(A T + Vxg)]^ 

^em = i(P L ) 2 + i[(n T +P|) 2 + (Vx Ax) 2 ] 

(47) 

with the first component representing the energy density 
of the matter field and the second component the energy 
density of the EM field. With these definitions the Hamil- 
tonian naturally represents the energy of the system as 
the sum of energies of the subsystems, rather than as the 
sum of superficially defined free and interaction energies. 

We note that the first term in J^em represents the 
energy density of the longitudinal EM field, while the 
second represents the energy density of the transverse 
EM field. The subsystems defined in this way are coupled 
because the velocity density "V and the electric field E 
do not commute; using for simplicity the Coulomb gauge 
(gT = 0), we obtain 

[r(x),^(x')] 
= [iV f (x)a • (V - ieA T (x))V(x) ) n TJ (x') + P Lj (x')] 

= ^ e? ^t( x ) a^ (\t (x _ x0 _ v . v ._±_) ^ (x) 

= ieift' (x)aiSij8(x — x')-0(x) = ijj(x)8(x — x') . (48) 

The (5-function ensures that ^(x) and E(x') are compat- 
ible observables for x/x'. Moreover it ensures that the 
matter field and the EM field energies are compatible in 
disjoint regions TZ 7 W C M 3 with K n W = 0: 

[J?S,#£m] = [ d 3 X [ dV[^4/(x),^EM(x')] =o. 

(49) 
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Using the commutator Eq. (|48|) and the commutator 
of the electric field and magnetic field energy 

i / "dV [B i (x),B i (x') 2 ] = ieyfeVj-B^x), (50) 

the equation of motion for the electric field is found to 
be 

E = VxB-j, (51) 

which is just one of Maxwell's equations. The remaining 
Maxwell equation 

B = -V x E (52) 

is found in a similar fashion as the equation of motion 
for the magnetic field. 

C. Energy and Causality 

In recent years a great deal of attention has been paid 
to the nature of energy transfer between separated ma- 
terial systems [20l - [25l |. specifically in the context of the 
two-atom problem of Fermi [l3j . 

In order to investigate causality at the microscopic 
level, Fermi considered two identical spatially separated 
atoms A and B. Initially atom A is energetically ex- 
cited while atom B is in its ground state and there are 
no photons present in the EM field. Einstein causality 
would appear to require that any changes in the energy of 
atom B be independent of atom A for all times less than 
the time it would take for a signal produced by atom A 
traveling at the maximal speed of propagation c to reach 
atom B. The most recent work concerning the Fermi 
problem was the proposal of a circuit QED experiment 
designed to test for any possible violations of causality 

The majority of theoretical proofs of causality in the 
Fermi problem involve using the bare states of a non- 
relativistic Hamiltonian in the Poincare gauge and the 
electric dipole approximation (EDA) [2l| - [24 |. The EDA 
dictates that the atoms in the Fermi problem couple to 
the Maxwell field at the c-number atomic center-of-mass 
positions. The problem itself can then be formulated in 
terms of a well defined, center-of-mass separation. More- 
over, the EDA ensures that the Poincare gauge dipole 
canonical momenta p, = mr 2 ; , i = A, B which define the 
bare atomic energies, arc purely kinetic [1, H|. 

At the same time in the Poincare gauge the field canon- 
ical momentum is identified as the (negative of the) lo- 
cal multipolar transverse displacement field 0, H| . Out- 
side the atoms, which in the EDA means away from the 
center-of-mass positions, this field coincides with the re- 
tarded electric field. The bare energies of the atoms are 
coupled through this field, which ensures that there are 
no violations of causality. 

A different proof presented in (2f| uses the abstract 
language of algebraic quantum field theory, relying quite 



generally on the primitively causal nature of relativistic 
quantum field theory resulting from the hyperbolicity of 
the relevant equations of motion. 

Here we will show using the gauge-invariant definition 
of M'm in Eq. (|47[) that changes in the energy density of 
the matter field at a point (t, x) are independent of the 
matter field at all points, which cannot be connected to 
(t, x) by a causal signal. The energy of the matter field 
in some closed region TZ C R 3 is merely 

ff5(t) = / d 3 xJP M (x,t). (53) 

Since 3%m is gauge-invariant the result does not rely on 
the use of a particular gauge and avoids any approxima- 
tions. 

We begin by calculating the equation of motion for 
J4?Mi which for simplicity is carried out in the Coulomb 
gauge 

jJ^m(x) = 

-\\ dV{[e^ t (x)a^ T) ,(x)^(x) ) n T , j (x / ) 2 ] 

+ [zV t (x)a l V^(x),P Lj (x') 2 ]} 

= y i (x)[S r , i (x)+£? L , i (x)] = y<(x)£((x). (54) 

Thus, together with the Maxwell equations (|51[) and (|52|) 
we have a system of equations, which can be written 

Jf M (*,t) = j(x,t)-E(x,i), 

j(x,t) = VxB(x,i)-E(x,t), 

□E(x,t) = -Vp(x,t)-j(x,t), 

□B(x,i) = V x j(x,t), (55) 

where □ is the d'Alembertian d 2 /dt 2 — V 2 . 

Now, first we note that using the second equation the 
first can be written in terms of the electric and mag- 
netic fields alone. Second we note that the remain- 
ing two equations are inhomogeneous wave equations for 
the (Cartesian components of the) electric and magnetic 
fields with source terms v(x,t) = V/o(x, t) +j(x, t) and 
fi(x.,i) = —V x j(x, t), respectively. These equations are 
hyperbolic and have well known retarded solutions of the 
form 

E(x,t) = E (x,t)+E r (x,t), 

B(x,t) = B (x,t)+B r (x,<) (56) 

where Eq and Bo satisfy the homogeneous equations 
□Eo = DBo = and E r and B r depend respectively 
on the sources is(x.'t r ) and £t(x', i r ) at the retarded time 
t T = t-\x- x'| [11, |3. 

Using these solutions, we can conclude that as desired 
Mm (x, t) at the point (t, x) depends on the matter field 
at points (t T ,x') only. 
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IV. CONCLUSIONS 

In this paper we have derived a Hamiltonian in an ar- 
bitrary noncovariant gauge, which could be taken as the 
starting point for the relativistic description of atoms and 
molecules. It also serves well as a means by which the 
gauge freedom of noncovariant canonical QED can be 
explored. We have discussed some implications of this 
gauge freedom and highlighted that a canonical partition- 
ing of the Hamiltonian is manifestly gauge-dependent. 

We have shown that a classical-type partitioning of the 



Hamiltonian in terms of velocities instead of canonical 
operators is gauge-invariant and suggested a possible ap- 
plication of such a partitioning in relation to the problem 
of causality in spatially separated material systems. 
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